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ABSTRACT 

We present six new transits of the exoplanet OGLE-TR-lllb observed with 
the Magellan Telescopes in Chile between April 2008 and March 2009. We com- 
bine these new transits with five previously published transit epochs for this 
planet between 2005 and 2006 to extend the analysis of transit timing vari- 
ations reported for this system. We derive a new planetary radius value of 
1.019 ± 0.026 which is interme diate to the previously reported radii o f 



1.067 ± 0.054 Rj (IWinn et all l2007f ) and 0.922 ± 0.057 Rj (IPiaz et al.l 120081 ). 



We also examine the transit timing variation and duration change claims of 



Diaz et al.l (l2008l ). Our analysis of all eleven transit epochs does not reveal any 
points with deviations larger than 2cr, and most points are well within la. Al- 
though the transit duration nominally decreases over the four year span of the 
data, systematic errors in the photometry can account for this result. Therefore, 
there is no compelling evidence for either a timing or a duration variation in this 
system. Numerical integrations place an upper limit of about 1 on the mass of 
a potential second planet in a 2:1 mean-motion resonance with OGLE-TR-lllb. 
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Introduction 



Transiting exoplanets provide a wealth of information for studies of the physical param- 
eters of planets and their environments. For example, the combination of several accurately 
timed transits of a known transiting exoplanet can be used not only to improve estimates 
of the planetary radius and orbital parameters of the star-planet system, but also to detect 
additional objects. Detecting potential variations of parameters such as the inclination and 
duration of the transits would indica te a precesing planetary orbit, potentially caused by 
another planet (IMiralda-Escudel 120021 ) . We can also use transit timin g to search for addi- 
tiona l planets or moons, as discussed in several recent theory papers ([Holman and Murray 



2ooa 



Agol et al. 



2ooa 



Hevl and Gladmaru 120071: iFord and Holmaru 12007 



Simon et al.l 



2007 



KippineJboOoi iKipping et al.l Eoool ) . The idea is that the presence of additional objects will 



perturb the orbit of the transiting planet, producing transit timing variations (TTVs) or 
transit duration variations (TDVs). Those TTVs and TDVs can be detected by monitoring 
transits over many orbital periods. The absence of such variations can be also used to place 
limits on the mass and orbital parameters of additional objects in those planetary systems 
and to gain insight into the systems' architectures. 

Recent observations show hints of timing variations for some transiting planets, but no 
definitive detection of additional planets or satellites has been reported using this technique. 
The most interesting results so far are (1) the absence of TTVs in several systems, which do 
not host planets more massive than several Earth masses in low-order resonant orbits (see a 
summary of constraints that can be pla ced in Table [1]); (2) th e tentative detection of TDVs 



in GJ436, roughly 3 minutes per year (jCoughlin et al 



20081 ). a trend consistent with the 



presence of a low-mass companion (< 12 Mq) in a close exterior b ut non- resonant orbit; thi s 
result is consistent with the 8 limit placed by transit timing (IBean and Seifahrtll2008l ): 
and (3) the preliminary detection of TTVs with a max imum residual of 156 ± 48 sec (3.3(t) 
over a period of 2 years reported by iDiaz et al.l (|2008[ ) for OGLE-TR-lllb, the subject of 
this paper. 

OGLE-TR-lllb is a 0.5Mj hot Jupiter orbiting its host star, a faint (/ = 15.5) K 
dw arf, every 4.01 days. This object was first announced as a transitin g planet candidat e 
by lUdalski et al.l (120021 ). and was confirmed to have planetary mass by jPont et al.l ( 120041 ). 
The physical parameters of the pl anet were refined over the next two years , with several 



new 



radial velocity measurements ( Gallardo et al. 2005; Silva and Cruzl 



20061 ) . The first high precision transit photometry was provided by IWinn et al. 



20061; 



two / -band transits of the planet on 2006 Feb 21 and Mar 5. Shortly after, iMinniti et al. 



Santos et al. 



(l2007h. with 



( I2OO7I ) published a F-band transit from 2005 April 9 an d noted that the rn idtime occurred 
5 minutes earlier than expected from the ephemeris in I Winn et al.l (120071 ) , although with 
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only t hree epochs they could draw no firm conclusions. A follow-up paper by iDiaz et al. 
( 120081 ) reported two consecutive /-band transits of OGLE-TR-lllb on 2006 Dec 19 and 23. 
Combining all five epochs, they concluded that the previously claimed TTVs were real, with 
the residuals spanning —156 ± 48 to +98 ± 39 seconds. Among other scenarios, they noted 
that if OGLE-TR-lllb were in an eccentric orbit with e ~ 0.3, the observed TTVs would 
be consistent with the pr e sence of an Earth-mass planet near an exterior 4:1 resonant orbit. 
Additionally, iDiaz et al.l (120081 ) noted two parameters with margi nally discrepant v alues 
across the fi ve tra nsits (see Table [2]). Compared to the results from lWinn et al.l (120071 ) . the 
Diaz et al.l ( 120081 ) values for the planetary radius disagreed at the 10% level, or 1.3cr, and 
the total transit duration differed by I.60". The radius ratio discrepancy was suggested to 
be the result of the parameters chosen for the image subtraction photometry, which focused 
on precise timing rather than on an accurate transit depth determination. The duration 
variation, if real, could be due to a perturber decreasing the orbital inclination, which would 
offer another way of determining the properties of the third body in the system suggested 
by their TTVs. 

Here we present six new transits observed during 2008 and 2009, which double the 
number of high-quality transit light curves available for OGLE-TR-lllb. In §[2] we describe 
the collection and analysis of the new data. In §|3]we describe the transit model fitting, and 
discuss additional sources of error not included in the formal fit. In §|l]we combine the six 
new transits with the five previously published observations and provide a new analysis of 
parameter variation in the OGLE-TR-111 system. In §|S]we discuss the implications of our 
results. 



2. Observations and Data Analysis 

All six new transits were observed between April 2008 and March 2009 in the Sloan i' 
filter with the new MagIC-e2v camera0 on Magellan. The MagIC-e2v camera has a field of 
view of 38" x 38" and a plate scale of 0.037" per pixel unbinned. With such high resolution 
and good average seeing at the site, blends are minimized and aperture photometry can 
be successfully applied even in fairly crowded fields. The camera can be operated in two 
different modes: single exposure mode, with a readout time of about 5 seconds per exposure, 
and frame transfer mode, with a readout time of only 3 milliseconds between frames in an 



^The MagIC-e2v detector, which shares a dewar with the older SiTe CCD, is identical to the red CCD 
on HIPO, a fast read-out direct imaging camera and o ne of the first generation instruments t o be flown 



on SO FIA; both cameras use the LOIS control software (jDunham et al.ll2004 iTavlor et aLll2004 lOsip et al 
20081 ). 
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image cube. Our first four transits were observed in single exposure mode. Tlie frame 
transfer mode first became available after engineering in July 2008, and was used for the 
last two transits of OGLE-TR-lllb. The gain and read noise of the first four transits were 
2.4 e-/ADU and 5.5 e- per pixel, respectively; after engineering, these values were changed 
to the current values of 0.54 e-/ADU and 5 e- per pixel. 

The exposure times during each transit were adjusted to maintain a minimum count 
level of about 10^ integrated photons, both for the target and multiple nearby comparison 
stars. For the 2008 transits, we collected unbinned (1x1) data with exposure times between 
30 and 120 sec, depending on the observing conditions, with an additional readout overhead 
per exposure of 5 seconds per frame. The 2009 data were collected in frame-transfer mode 
with the camera binned 2x2, which yielded an improved sampling rate of 15-30 sec per frame. 
Details of the observing settings are noted in Table Ei 

Accurate timing is of the utmost importance for this project, so special care was taken to 
ensure that the correct times were recorded in the image headers. For the 2008 transits, the 
start times for each image were recorded from a network time server, which was verified by 
eye to be synchronized with the observatory's GPS clocks at the beginning of each night. For 
the 2009 observations, the times came from a PC104 (a small embedded control computer), 
which received unlabeled GPS pulses every second. As with the network time server, the 
PC104 was synchronized with the observatory's GPS before each transit observation. In 
both cases the time signals written to the image headers agree within one second with the 
GPS time. During both 2009 transits, a software failure caused the times for a few image 
cubes to not be recorded in the headers of the images, but we were able to reconstruct 
the observation times with precisions better than a second from detailed system logs. One 
second is a conservative estimate of the intrinsic error for the start time for each frame, and 
is significantly smaller than the mid-transit times errors. 



2.1. Data Analysis 

All data were overscan corrected and flattened using IRAF@ The photometry was per- 
formed using the IRAF routine phot, part of the apphot package. Depending on the binning 
applied and the seeing during each transit, a subset of apertures between 6-25 pixels in ra- 
dius were examined for the target star and each comparison star. Between one and seven 



^IRAF is distributed by the National Optical Astronomy Observatories, which are operated by the As- 
sociation of Universities for Research in Astronomy, Inc., under cooperative agreement with the National 
Science Foundation. 
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stars were selected as comparisons from the list of 10-20 field stars present in each frame. 
The comparison stars had to pass several selection criteria: to be similar in brightness to 
the target, to not be blended, and to not be variable. The best apertures were identified as 
the ones which yielded the smallest scatter in the out-of-transit flux. The choice of position 
of the sky background annulus also influenced which stars could be used for comparison; see 
§ 13.31 for a discussion of systematic errors resulting from aperture settings. We explored sky 
regions with inner radii from 20-40 pixels and 10-30 pixels width, and selected the one which 
(a) provided several comparison stars of suitable brightness that were not variable, and (b) 
produced the lowest noise in the out-of-transit baseline while not introducing spikes or other 
obvious problems in the transit light curve. 

In each of the resultant light curves, the out-of-transit baseline was examined for linear 
correlations with several variables: airmass, seeing, telescope azimuth, (x,y) pixel location, 
and time since beginning of transit. The parameters chosen are either directly correlated 
with physical phenomena that affect photon rates (e.g. seeing, airmass), or are proxies 
for other effects (e.g. the telescope azimuth is correlated with the de-rotator rates, which 
were not recorded for several transits). Only the 2009 transits exhibited significant trends, 
against telescope azimuth (both) and seeing (20090313). For each transit, we corrected 
those trends by successively subtracting linear fits to each variable. Detrending was critical 
in producing usable light curves for these two transits, but may also have introduced smaller- 
order systematic effects, particularly in the transit depth (see also § 13. 3p . 

We now briefly describe the observations and the photometric reduction of each transit 
dataset. All transit fluxes and times are available online; an excerpt is shown in Table |H 

2.1.1. 2OO8O4I8 

Transit 20080418 was observed during engineering time just after the e2v CCD was first 
installed on the telescope; due to engineering constraints, only the second half of the transit 
was observed. The field was repositioned before egress to eliminate diffraction spikes from 
a nearby bright star by moving it further off-chip. The airmass was low and fairly constant 
(1.2-1.3) and the seeing was good, increasing slightly from 0.4" to 0.5" during the transit. 
The out-of-transit data showed no apparent trends with the parameters checked. 
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2.1.2. 20080422 

Transit 20080422 was observed on an intermittently cloudy night with highly variable 
transparency, with counts on the target star varying by a factor of 6 within a few frames. 
We found that eliminating the lowest count frames, those with fewer than 300,000 counts 
on the target star, significantly decreased the scatter of the light curve. The seeing ranged 
from 0.5-0.6", and the airmass was low and fairly constant (1.2-1.3). The out-of-transit data 
showed no apparent trends with the parameters checked. 

2.1.3. 20080512 

Transit 20080512 had stable photometric conditions for the entire pre-transit baseline. 
During the transit there were two drops in target counts (by a factor of 2) that coincided 
with sudden seeing jumps (0.4" spiking to 0.6"). The field also drifted substantially (by 
~100 rows and ~100 columns) due to tracking problems; about 30 minutes of post-transit 
baseline had to be discarded because of strong image elongation. The airmass ranged from 
1.2 to 1.7. The out-of-transit data showed no apparent trends with the parameters checked. 

2.1.4. 20080516 

Transit 20080516 had very stable photometric conditions for most of the transit. The 
seeing gradually increased from 0.4" to 0.6", and the airmass ranged from 1.2 to 2.0. The star 
also drifted substantially toward the end of transit (by ~ 200 rows and ~ 200 columns) for 
unknown reasons. The out-of-transit data showed no apparent trends with the parameters 
checked. 



2.1.5. 20090217 

Transit 20090217 was the first of OGLE-TR-lllb to be observed with the new frame 
transfer mode. The seeing fluctuated from 0.7-1.1", while the airmass decreased from 1.8 to 
1.2. This transit was detrended for a slope correlated with the telescope azimuth. 
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2.1.6. 20090313 

Observations for transit 20090313 began late due to telescope problems and thus there is 
no pre-transit baseline. Seeing conditions were initially very poor, spiking to 2", but improved 
substantially during the second half of the transit, to around 0.6". The airmass ranged from 
1.2 to 1.5. The huge disparity in seeing made photometry on this transit challenging, and 
the best light curve resulted from using a small aperture around the target and only the 
brightest single comparison star, which saturated for a few frames that had to be discarded. 
This transit was detrended for slopes correlated with the seeing and the telescope azimuth. 



2.1.7. Literature light curves 

To eliminate any uncertainty in comparing our transits to previously published transits, 
which might result from different models, fitting procedures, or fixed parameter values, we 
have obtained tables of the times and fluxes for each of the five transits drawn from the 



literature: 20050409 flMinniti et al. 20071) . 20060221 and 20060305 flWinn et al.l 120071 1. and 
20061219 and 20061223 (jPiaz et al.ll2008l ). We have used the original photometry, except for 
converting from magnitudes into fluxes (where applicable) and converting the mid-exposure 
times from Heliocentric Julian Day, or HJD, into Barycentric Julian Day, or BJD, to be 
consistent with the rest of our analysis; the difference between the two time systems is much 
smaller than our errors, though, a few seconds at most. All values presented for the literature 
light curves are taken from our re-analysis of the published photometry using our model and 
fitting, unless otherwise noted. 

The light curves for all six new transits and the five literature transits are shown in 
Figure [H together with the best joint model fit, which will be described in detail in § |3l 



3. Transit fitting results 



3.1. Model 



Each light curve was fitted with the iMandel and Agoll (|2002[ ) algorith m to generate an- 
alytic al models, using the basic optimized model-fitting code described in ICarter and Winn 
( I2OO9I ). but without the wavelet analysis. In the models we assumed that OGLE-TR-lllb 
has zero obliquity, oblateness and orbital eccentricity. We used a quadratic limb darkening 
law of the form 

/(r) = 1 - Mi(l - Vl - r2) - U2{1 - Vl - r2)2, (1) 
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with the in itial parameters for ui and U2 set to the values for the appropriate filter (j Claret 



2OOOI . l2004j ). Although analyses of differen t limb darkening laws have shown that using a 



non-linear limb darkening law is important (ISouthworthll2008[ ) . it is generally not possible to 
fit both quadratic limb-darkening coefficients except on the highest quality, typically space- 
based, data. We thus fixed the quadratic term U2 and only fit for the linear term, ui. (We also 
fixed Ml for the sparsely-sampled transit on 20050409, the onl y one observed in V band.) The 
values for ui and U2 are calculated using the jktld program by Southworth (2008]§, assuming 
T = 5044 K, logg = 4.25 cm/s^ [M/H] = 0, and Vrmcro = 2 km/s. We used the limb 
darkening v alues corresponding t o the Sloan i' filter for the new data, t he Vg filter f or the 
tra nsit from iMinniti et al.l ( 120071 ). and the Ic filter for the transits from lWinn et al.l (|2007| ) 
and Diaz et 817 ( 2008 ). as listed in Table |5] (though we note for completeness that the actual 
filters used in the literature light curves. Mould V, CTIO / and Bessel /, respectively, do not 
correspond precisely to the Cousins V and / filter parameters that were available). We fixed 
the orbital period to P = 4.01445 days; later experiments with slightly different values had 
little effect. The other free parameters in the model are the radius ratio, k, inclination, i, 
semimajor axis in stellar radii, a/R^,, out-of-transit flux, Fqot and transit midtime, Tg- We 
assume throughout the fits t hat = 0-§l ^* = 0-^3 -R©, and Mp = 0.52 Mj, based 
on the spectroscopic work of ISantos et al.l (l2006l ). 



3.2. Light curve fits 



To determine the best fit value and error of each model parameter , we u sed a Monte 
Carlo Markov Chain (MCMC) method, as described in lCarter and WinnI (120091 ). The initial 
values for each parameter were computed by a joint least squares fit to each light curve 
independently. We then weighted each light curve by the reduced of fhis fit so that the 
new reduced = 1- (By doing this we are assuming the transit model is correct in order to 
determine the error on each transit, rather than assuming a noise model, e.g. photon noise, 
in order to test the transit model.) A joint least squares fit of all weighted transits had a 
reduced = 1-05. Starting from the initial least squares values, we constructed chains of 
1,000,000 links, where the acceptance rate for each parameter is between 20-60%. We fit 
all eleven transits simultaneously, assuming common values for k, i, a/R^,, and ui^x (where 
X is the appropriate filter), with U2,x fixed, but fitting each transit for its own Fqot and 
Tc- We did not fit for an airmass slope (see § 13.31 The first 50,000 points of the MCMC 
were discarded to eliminate bias toward the initial conditio ns. We created three inde pendent 
MCMC chains, checked that the Gelman-Rubin statistic (jCelman and Rubinlll992l ) is close 



•^http: / /www. astro. keele.ac.uk/ jkt/ codes /jktld. html 
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to 1 to ensure convergence, and then combined the chains to determine the distribution of 
all parameters, including the total duration of the transit, T14 (the time from first to fourth, 
or final, contact), and the impact parameter, b, which are derived from fitted parameters. 
We plot the best model fit with the data in Figure [T] and tabulate the fit results in Table El 
where we report for each parameter the median value and the 68.3% credible interval (the 
equivalent to a la standard deviation if the distribution is Gaussian). The distributions for 
each parameter are shown in Figure [2l 

The new radius ratio for OGLE-TR-lllb based on an analysis of all eleven light curves 
yields a planetary radius Rp = 1.019 ± 0.007 Rj, if we consider only the formal fit errors; 
accounting for the error on the stellar radius, which is now the dominant source of error, we 
find a more realistic error bar is Rp = 1.019 ±0.026 Rj. Note that if we only use the six new 
light curves, which have more consistent radii, the value for the radius ratio is very similar 
(formal fit Rp = 1.015 ± 0.009 Rj, or Rp = 1.015 ± 0.026 Rj with stellar errors). 

As a test of the robustness of our parameter determination, we also ran additional 
MCMC fits for each transit independently, with results in Table O Although most of the 
parameters agree within the formal la errors between the individual and joint fit s, there are 



some notable exceptions. In § 14.21 we investigate variability (previously noted by iDiaz et al. 



( 120081 )) in both the radius ratio, k, and the total transit duration, T14. 



3.3. Systematic Errors and Correlated Noise 

One of the most apparent results from the fits to individual light curves was that the 
radius ratios are similar for transits observed on the same instrument and reduced by the 
same group. This may indicate a degree of subjectivity in the light curve generation process, 
both from the choice of photometry method (e.g. aperture, image subtraction, deconvolution, 
etc.) and from the specific choice of reduction parameters (e.g. aperture size and sky region 
for aperture photometry). These choices can result in systematic errors in the transit depth, 
particularly when comparing transits from multiple sources. Both image subtraction and 
aperture photometry require fine-tuning a number of parameters, and there is no single 
prescription for how to get the absolute best light curve: the same method applied to the same 
transit could produce similar quahty hght curves, as measured by the scatter of residuals or 
out-of-transit flux, whi ch nonethele s s diff er in depth by more than the formal fitted errors. 



It has been noted by IWinn et al.l (120071 ) that with image subtraction, slight changes in 
both the difference flux and the reference flux can cause the measured radius ratio to vary 
by a few percent of its value, although their estimate of that effect on their own data, 
6k = 0.0002, is much less than our formal fit error of 0.0008. This effect was also alluded 
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to by iDiaz et al.l ( 120081 ) as an explanation for their shallow depths compared to previous 



results, althoug h they did not prov ide a numerical estimate of the magnitude of this effect. 



An analysis by iGillon et al.l ( 120071 ) of a different transiting planet, OGLE-TR-132b, found 
that image subtraction is particularly prone to misestimating the transit depth, compared to 
the alternative methods of aperture and deconvolution photometry. This effect for OGLE- 
TR-132b in their data causes the radius ratio to differ by 1-2% depending on the choice of 
parameters . If si milar levels of error were present for OGLE-TR-lllb, particularly in the 



Dfaz et al.l (120081 ) curves which were acknowledged to not be optimized toward finding the 
correct depths, the systematic error on the radius ratio would be 0.0013 — 0.0025, comparable 
to the formal fit error of 0.002 on the individual curves. Assuming a median value of 0.002 
for the systematic error and adding it in quadrature with the formal fit error, a better error 
estimate on the radius ratio would be 0.003. 

We attempted to quantify the systematic error for aperture light curve generation as 
follows. For each of our light curves, we used 4-5 sets of apertures, sky radii and widths, 
and different comparison stars, with a goal toward minimizing the scatter in the out-of- 
transit baseline. The choice of comparison stars in particular depends on the aperture and 
sky choices, because certain stars are usable under some choices but not others, especially 
stars of very different brightnesses. Additionally, we examined the effect of detrending light 
curves, which must be done carefully because systematic trends in the data can distort the 
measured radius, but so can an incorrectly-removed slope. We only detrended transits which 
had strong slopes in the out of transit baseline or otherwise had distorted shapes, and then 
only for the parameters which best corrected the shape defects. It is possible that residual 
trends against parameters we did not consider, or slightly nonlinear trends, could remain 
in the data and distort our estimate of the radius; we chose to stick with linear trends 
against a few meaningful physical parameters to avoid introducing unnecessary complexity. 
To test whether any important slopes remained, we ran a joint MCMC fit which includes a 
differential extinction term (a trend with airmass), but found that the fitted airmass slopes 
were slight and the difference in resulting parameters were in all cases less than the formal 
la errors. The results of our explorations of both the aperture and sky choices and the 
detrending parameters are that we can produce light curves with similar shapes and scatter, 
and that the radius ratios vary by 0.001-0.004 for an individual transit. (Some transits are 
much more resilient to parameter choices than others.) Thus, for aperture photometry also 
the systematic error on the radius ratio based solely on the parameter choice is of order the 
formal fit errors. 

To estimate the amount of correlated noise in the light curves, we used two methods. 
The first is residual permutation, which shifts the residuals for each transit through every 
point in time and adds it to the best model fit; we also assumed time invariance and reversed 
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the residuals, then permuted again, for a total of 214-1600 curves for each transit. We fit 
a least-squares transit model to each permuted curve for all eleven transits. For the radius 
ratio and the transit midtime, the errors from the residual permutation method for both 
values were greater by a factor of 1-3, depending on the light curve (note that this factor is 
not the same for each parameter; see Table [6]). We got similar results when we ran a joint 
fit of all 11 transits with 10,000 curve ensembles, randomly selecting for each transit one of 
its individually permuted light curves. 

An alternative way to estimate the erro r contributed by correlated noise is to calculate 



how the noise scales with time averaging (jPont et al.l l2006l ). We calculate the standard 



deviation on the residuals in bin sizes from 10-30 minutes and compare that value to what 
we would expect if the noise behaved like Poisson noise (i.e., a decrease in the noise with 
a/ZV points). We calculate the amount by which the real noise is greater than the estimated 
noise, and find that it is greater by a factor of 1.5-3 times the purely Poisson noise level, 
depending on the transit. The increased noise factors agree with the values found by residual 
permutation, and for simplicity we use the scaled errors from the residual permutation 
throughout. 

For all of these reasons, the formal fit errors reported in Table |5] are underestimated due 
to both systematic errors and correlated noise. The errors we adopted have been inflated 
based on the residual permutation method. This means that, once systematics are included, 
neither the observed timing variations for 20050409 and 20080418 nor the duration variation 
are statistically significant (see § 14.11 and 



Results 



4.1. Timing 



The central midtimes for all 11 transits that we fit a re summarized in Table [5] and 
illustrated in Figure IH Recently iPietrukowicz et al.l (120 lOf ) have reanalyzed the photome- 
try for 30 transits of OGLE planets and planet candidates, among t hem OGLE-TR-lllb , 
and they have found a different midtime than originally reported in iMinniti et al.l ( 120071 ): 
Tc,new = 2453470.5676 ± 0.0005, compared to the original published value of Tc\orig = 
2453470.56413 ± 0.00067, a difference of 300 seconds. (It is not clear what is the source 
of such a large shift, but one possibility is a mistake in the UTC-BJD correction.) Signifi- 
cantly, the new time is much closer to the expected time of transit. 

Another potential pitfall when comparing times from multiple groups has recently noted 
by Eastman and Agol (in prep). Most researchers, and indeed most common conversion tools 
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(e.g. barycen.pro in IDL and setjd in IRAF) by default omit the correction from UTC to 
TT, which in 2009 was 6 6.184 seconds. We hav e confirmed that the times pubhshed by 
Pietrukowicz et aLl (120101 ) and IWinn et al.l ( 120071 ) do not account for the UTC-TT co rrec- 
tion (personal co mmunications), and we assume that the times in iMinniti et al.l (120071 ) and 
Diaz et al.l ( 120081 ) likely did not either. We have therefore added the appropriate correction 
to the reported BJD times for these light curves. (Note that the smaller order deviations 
introduced by using UTC rather than TT times in calculating the BJD correction terms are 
at most a few seconds, and for this work those deviations fall well within the timing errors; 
however, with higher precision data on other systems it would be very important to consis- 
tently calculate the BJD times.) All of the transit midtimes in Table [5] have been corrected 



to the BJD-TT system; additionally, we have added 64.184 seconds to the lPietrukowicz et al. 
(I2OIOI ) midtime to get Tc = 245470.56834 (BJD), the value used in all subsequent calcula- 
tion. 



The top panel of Figure H] shows the ephemeris from iDiaz et al.l (120081 ) , derived from 
the first five transits, along with the timing residuals for all eleven transit based on our 
joint-fit values (Tabled]). To correct for the linear drift in the residuals, we calculate a new 
constant-period transit ep hemeris, omitting the half -transit (20080418), which has very large 
errors, and using the new iPietrukowicz et al.l (120101 ) time for 20050409, and find: 



Tc = 2454092.80717(16) [BJD] + 4.0144463(10) A^. (2) 

where Tc is the predicted central time of a transit, and A^ is the number of periods since 
the reference midtime, and the values in parentheses are errors on the last digits. We find 
almost identical ephemeris values if we use all transits and if we use the original time for 
20050409, although the errors are several times greater in both cases. Our adopted fit has a 
reduced = 0.5. 

The lower panel in Figure H] shows the new ephemeris and the timing residuals. (Note 
that the same times are used in both panels, and only the ephemeris has changed.) With la 
errors ranging from 36 s to 114 s, only the original 20050409 time is more than 2a from zero, 
and of the other transits only the half transit 20080418, which is inherently less t rustworthy, 



i s mo re than la. Thus, we conclude that the timing deviations reported by iDiaz et al. 



(I2OO8I ). which depended heavily on the old time for 20050409, do not exist, and we see no 



evidence for timing variations in our data. 
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4.2. Parameter variation 



Diaz et al.l (|2008() found that thei r value for the total duration was 4.4 minutes shorter 
than that found by I Winn et al.l (120071 ). a 1.6a result given the respective quoted errors. If 
this decrease is real, it would be of great interest, since a likely explanation would be that the 
inclination of OGLE-TR-lllb is precessing, possibly due to the presence of another planet. 
On the other hand, the variation could be due to errors in the photometry or undetected 
correlated noise. Our values for the best fit duration for each transit are plotted in Figure [3] 
and tabulated in Table [6l along with the parameters the duration was derived from: i, k, 
and a/R^. Note that the errors in this table have been increased from the formal fit errors 
by a factor derived from the residual permutation method discussed in § 13.31 

At first glance there does appear to be a decrease in duration over time. To compare 
the durations of the 10 full transits (excluding 20080418) we fit the data using two models: 
(1) a flat line, corresponding to a constant duration, with reduced = 0-9, and a value of 
9636 ± 58 s, very similar to the joint-transit fit value, and (2) a sloped line, with reduced 
= 0.6, and a slope of —0.24 ± 0.12 seconds per day. This is a much shallower slope than 
the hint of a trend reported by iDiaz et al.l (120081 ). and is not significant given the good fit 
achieved with the fiat line. 

Because the duration is a derived quantity, we must examine the parameters on which it 
depends {i, a/R^, and k). Both the inclination and the semi-major axis exhibit slight slopes, 
but with low significance given the errors {i changes by 0.22 ± 0.09 degrees per year with 
reduced = 0.2, and a/R^ by 0.0007 ±0.0002 stellar radii per year with reduced = 0.2). 
If the duration were really decreasing because the planet is precessing, this would be due to 
the planet moving away from the center of the stellar disk and hence the inclination would 
decrease; instead, the best-fit slope for the inclination is slightly positive, another indication 
that we are not picking up on a real effect. 

On the other hand, the radius ratio does have real variations between transits, given the 
curre nt photometry. W e find variation from a low of A; = 0.118 ± 0.002 (transit 20061219, 
frombiaz et al.ll2008h to a high of k = 0.132 ± 0.002 (20060221 and 20060305, both from 



Winn et al 



20071 ). with the rest of the transits in between ( see Table [6]) . (No te that our 



individual fits for the radius and error of the two transits from Winn et al 



(120071 ) agree with 



those cited in that paper, indicating that our fitting methods are comparable.) The best fit 
line with a slope is not statistically significant (within la of 0), and the most likely expla- 
nation of the variation in radius depths is due to systematic effects in how the photometry 
was reduced (see § 13.31) . 



If the star were active, the presence of stellar spots or active regions can affect the 
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observed ra dius depths ( s ee e.g . IPont et all (120081) for HP 189733. iRabus et al.l f l2009h for 
TrES-1, and lCzesla et"al] toQ^ : Ifluber et al.l feoooh : Isilva-Valio et al.l koi(h for CoRoT-2, 
all of which are known or theorized to be spotty stars). There is no record of variability 
for this star in the current literature. We examined two publishe d data sets with observa- 



tions over multiple non-transit nights: 4 nights of VIMOS data (jPietrukowicz et al. 



Minniti et al.l 120071 ) and data spanning 115 nights from the OGLE survey 



2010 



Udalski et al. 



( I2OO2I ). We found that the long-term flux was stable to within a few mmag in both datasets 
(3 mmag and 5 mmag respectively). 

Since the trends in the parameters i, a/R^,, and k are slight and unlikely to be physical, 
we cannot conclude based on the available data that the observed duration variation is a 
real effect. 



4.3. Limits on perturber mass 



Although there is no clear evidence of TTVs beyond the 2a level in the current dataset, 
we can use the TTVs reported in Table [7] and shown in Figure H] to place upper limits on 
the mass and orbital separation of a hypothetical perturbing plane t in the system. For tha t 
purpose we use an implementation of the algorithm presented in ISteffen and Agoll ( 120051 ). 
kindly provided by D. Fabrycky. 

We explored the full perturber's mass parameter space for interior orbits and exterior 
exterior orbits for orbital periods from 0.9-17.5 days (0.2-4.4 times the orbital period of 
OGLE-TR-lllb), and small initial eccentricity Cc = 0.05 (we also examined Cc = and 
Cc = 0.3). All the orbits were a ssumed coplanar and the orbi tal instability regions in each 
case were determined following iBarnes and Greenberd (120061 ). The orbital period of the 
perturber was increased by a factor of 1.1 for each step. For each period, the mass of the 
perturber also increased from an initial value of 0.0001 until reaching a mass that would 
produce a TTV equivalen t to the 3a confide r ice le vel of our results (i.e. Ax^ = 9). We 
used the time reported by iPietrukowicz et al. (I2OIQI) for the transit 20050409 instead of the 
value we fit from the original photometry from iMinniti et al.l (120071 ) , since the revised value 
is more in line with expectations; simulations run omitting that transit yield similar results. 
The mass limits placed by these tests are illustrated in Figure HI 

The constraints placed on the perturber's mass are strongest near the low-order mean 
motion resonances, particularly in the interior and exterior 2:1 resonances, where we are 
sensitive to objects as small as 1 M® and 0.6 M®, respectively if Cc = 0.05 (cc = yields 
even smaller constraints, but we choose to cite the more conservative value). No meaningful 
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constraint c an be placed on the region of the 4:1 external resonance, which was identified by 



Diaz et al.l (|2008[ ) as a possible location for a 1 perturber that could explain their TTV, 



but since we do not reproduce their TTV such a perturber is no longer necessary. 

Finally, our 0-C data, with a la precision of 36-114 seconds (after accounting for system- 
atic errors), cannot constrain the presence of moons aro und OGLE-TR-lllb , which would 
introduce a TTV signal of 4.6 seconds for a 1 moon ( iKipping et al.l 120091 ). 



5. Conclusions 



We have tested the previously claimed presence of TTVs and TDVs in the OGLE-TR- 
111 system by adding six new transit epochs, obs e rved between 2 008 a n d 2009, to t he fiv e 
previously published results bv Iwinn et al.l (120071 ): Minniti et al.l fcoOTH : IPiaz et al.l fcoosl l. 
This new analysis not only doubles the number of available data points, but also extends 
the TTV baseline from two to four years. In addition, combining the six new transits data 
allows us to provide a new, more precise value of the radius of this planet. We find a new 
radius f or the planet of 1.0 19 ± 0.026 -Rj, which is intermediate to the previously reported 
radii by lWinn et al.l (120071 ) and iDfaz et al.l (|2008[ ). and is more precise. 



We find a slight variation over time of the duration of the transits of OGLE-TR-lllb, 
as well as variations of other parameters, such as the inclination and semimajor axis of the 
orbit. Those variations could, in principle, be attributed to perturbations of the orbit of 
OGLE-TR-lllb produced through interaction with additional planet (s) in the system, but 
we demonstrate that the variations can be instead explained by systematic errors in the 
data, and therefore should not be attributed to other planets. 

We have also computed the transit midtimes of our new transits with formal precisions 
of 20-40 seconds, and more accurate precisions of 35-50 seconds for the full transits (and 
almost 2 minutes for the half-transit) once systematic errors are considered. The errors 
on the literature transits similarly increased when the photometry is refit using the same 
method to account for systematics, to 60-110 seconds depending on the light curve. 

A longer time baseline and more precise timing data is still necessary to test further for 
the presence of other planets in the OGLE-TR-111 system, especially in potentially stable 
non-resonant orbits, but with the present results we conclude that OGLE-TR-111 belongs 
in the category of systems summarized in Table |T] for which there is no sign of additional 
planets more massive than a few in low-order resonant orbits, including a limit of 1 
near the 2:1 resonances. The presence of massive (Earth-like) moons around OGLE-TR-lllb 
is still possible, but to detect those we would require timing precision of a few seconds or 
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better, beyond the current capability of ground-based instrumentation for this system. 
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Fig. 1. — Eleven transits of OGLE-TR-lllb. All available high-quality light curves are 
plotted vs. orbital phase, with the data binned to 2 minutes to aid comparison. The joint 
model fit (solid lines) were calculated using the parameter values in Table [Sj the stated 
standard deviation is the residuals from the joint model fit. Table H] shows the unbinned 
data; a full table is provided online. 
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Fig. 2. — Smoothed histogram of normahzed parameter distributions, from which the pa- 
rameters in Table \5\ are derived. The sohd hne is the median value (which is very close to 
the mean value in all cases). Note that several of the distributions, particularly a/i?^, and 
6, are not strictly Gaussian. The dashed lines show the 68.3% credible interval. These values 
were calculated for 2,850,000 links. 
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Fig. 3. — Parameter variation vs. time for all transits of OGLE-TR-lllb, based on individual 
MCMC fits (Table [6]), for (clockwise, from top-left): k, i, a/R^,, and T14. The errors have 
been scaled upward based on the factor calculated from residual permutation. The values 
derived from the joint MCMC fit to all transits (Table [5]) are plotted as solid black lines with 
±lcr errors. The dashed red lines indicate the best sloped fit with ±lcr errors, although all 
fits are only marginally significant (within la of a constant value for the radius ratio and 
within 2cr for the other parameters). 
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Fig. 4. — Observed minus calculated midtimes for O GLE-TR-ll l b. T op panel: timing 



residuals for eleven transits using the ephemeris from iDiaz et al.l (120081 ). Bottom panel: 
timing residuals using the new ephemeris (Equation |2j The 0-C values and errors are 
shown in Table U\ and were calculated using the formal fit midtimes reported in Table [5] 
but rescaling the errors to more realistically account for systematic noise (see § 13. 3p . The 
solid line represents zero deviation from expected time of transit, while the dashed lines 
represent the la and 3a errors on the calculated orbital period, indicating the slopes that 
result for a mis-determi ned period. We plot our calculated midtime for 20050409, based 
on the photometry fr om iMinniti et al. J2007h. as an open circle, and the revised midtime 
and error reported by iPietrukowicz et al.l ( l2010h as a s olid triangle. The ri e w eph emeris was 
calculated using only the solid symbols, i.e., using the IPietrukowicz et al.l (l2010l ) time. 
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Fig. 5. — Upper mass limit on potential companion planets vs. period of the perturber. We 
examine companions with initial Cc = 0.05(black) and Cc = 0.0 (gray). The constraints are 
strongest near the 2:1 mean-motion resonance, where objects as small as 1 would have 
been detectable; other interior and exterior resonances are also labeled. The shaded grey 
region shows the instability region for a 1 Mq, following iBarnes and GreenbergI ([20061). 
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Table 1. Maximum Allowed Perturber Mass (M^) Constrained by Transit Timing 



System 


Max O-C 


1:2'' 


2:3 


Interior 


3:2 2:1 


Exterior 


Reference 




(sec) 






(non-resonant) 




(non resonant) 





CoRoT-1 


78 ±24 


unstable 


unstable 


unstable 




4 


200 (P < 4d) 




< 60 










1 


10 (P < 3.5d), 
















100 (P < 5d) 


GJ 436 


< 60 












8 (P < 5d) 




< 120 


10 


2 




2 


20 


20 (P < 5d) 


HAT-P-3 


< 60 


0.7-1.0 








30-40 




HD 189733 


< 45 


4 


1 




8 


20 


32 (P < 5d) 


HD 209458 


70 ± 50 


0.3 




20 


0.3 


0.3 


















17 (P < 7d), 
















100 (P < lOd) 


TrES-1 


107 ±17 


1 


2 


3 


2 


1 


100 (P < 6d) 


TrES-2 


257 ± 27 


2 


15 


30 




1 


50 (P < 7d) 


TrES-3 


70 ±30 


3-4 








10-15 





Bead ll2009f) 



Csizmadia et al.l ||2009|) 



Bean and Seifahrj |2008h 
Ballard et all ll2009h" 
Gibson et aU ll20q9bl'l 



Miller^Ricc^_et_^lJ 

Miller-Ricci et al.l 



2008^ ) 
2008aii 



Agol and SteffenI | |2007^ 



Steffen and 


AgoJ 




Rabus et al. 




200£ 


) 


Rabus et al. 




200£ 


) 


Gibson et al.l ll2009al) 



iooi); 



''A perturber in an n:m resonance completes m orbits while the known planet completes n. 



Table 2. Selected Parameters for OGLE-TR-lllb 



Radius" 

(Rj) 


Midtime 
(JD) 




Duration 
(hours) 


Impact parameter 


Filter 


Reference 


0.97 ±0.06 


24552330 


44867 (fixed) 




- 0.68 


/ 


Pont et al. (2004); Santos ct al. (2006) 


1.01 ±0.06 


24553470 


56389 ± 0.00055 


2.9 ±0.1 


- 0.65 


V 


Minniti ct al. (2007): Di'az ct al. (2008) 


1.067 ± 0.054 


24553787 


70854 ± 0.00035 


2.743 ±0.033 


0.25 - 0.55 


I 


Winn et al. (2007) 


1.067 ± 0.054 


24553799 


75138 ± 0.00032 


2.743 ±0.033 


0.25 - 0.55 


I 


Winn et al. (2007) 


0.922 ± 0.057'= 


24554088 


79145 ± 0.00045 


2.67 ±0.014 


0.38 ±0.2 


I 


Di'az et al. (2008) 


0.922 ± 0.057*= 


24554092 


80493 ± 0.00045 


2.67 ±0.014 


0.38 ±0.2 


I 


Di'az et al. (2008) 


1.025 ± 0.047 


24553470 


56385 ± 0.00100 


2.75 ±0.11 


0.39 ±0.28 


V 


this work"* 


1.066 ± 0.038 


24553787 


70860 ± 0.00076 


2.73 ±0.047 


0.30 ±0.14 


I 


this work"^ 


1.062 ± 0.051 


24553799 


75135 ± 0.00090 


2.76 ±0.065 


0.39 ±0.16 


I 


this work"^ 


0.952 ± 0.076 


24554088 


79118 ±0.00138 


2.62 ±0.060 


0.38 ±0.12 


I 


this work"^ 


0.968 ± 0.064 


24554092 


80494 ± 0.00129 


2.66 ±0.071 


0.35 ±0.15 


I 


this work"* 


1.045 ± 0.036 


24554574 


54805 ± 0.00078 


2.40±0.393<= 


0.51 ±0.27 




this work"* 


1.030 ± 0.029 


24554578 


55395 ± 0.00061 


2.73 ±0.062 


0.36 ±0.19 




this work"* 


0.996 ± 0.022 


24554598 


62680 ± 0.00047 


2.63 ±0.037 


0.22 ±0.11 




this work"^ 


1.028 ±0.021 


24554602 


64098 ± 0.00053 


2.69 ±0.051 


0.30 ±0.12 




this work"^ 


1.006 ± 0.025 


24554879 


63787 ± 0.00056 


2.68 ±0.054 


0.28 ±0.12 




this work"* 


1.025 ± 0.022 


24554903 


72515 ± 0.00039 


2.61 ±0.058 


0.16 ±0.11 


i' 


this work"* 



''Except as noted, assuming P* = O.83P0 

•^Using P, = O.8IIP0; if P* = 0.83Po is used, the radius is Rp = 0.944Pj 
■^Half-transit 

■^AU values from this work are taken from individual MCMC fits to each light curve, with with errors infiated to account for 
correlated noise; see 13.31 
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Table 3. Observational and Aperture Photometry Parameters for Six New Transits 



Transit Frames Exp. Time Binning Readout Ncomp Aperture* Sky radius Sky width Scatter'' Est. Poisson 



(UT) 


used (discarded) 


(sec) 




(sec) 




(pixels) 


(pixels) 


(pixels) 


(mmag) 


(mmag) 


20080418 


178 (0) 


30-60 


1x1 


5 


3 


12.8 


35 


15 


1.3 


1.1 


20080422 


110 (21= ) 


120 


1x1 


5 


4 


16.4 


30 


10 


2.0 


0.8 


20080512 


241 (20<* ) 


60 


1x1 


5 


3 


19.2 


30 


15 


1.5 


0.8 


20080516 


276 (4« ) 


30-100 


1x1 


5 


7 


17.4 


30 


10 


1.5 


0.8 


20090217 


800 (0) 


30 


2x2 


0.003* 


6 


19.2 


25 


20 


1.2 


1.1 


20090313 


600 (428 ) 


15-30 


2x2 


0.003^" 


1 


9.6 


20 


30 


1.5 


1.5 



^Radius around star. 

''Standard deviation of the residuals on data binned to 120 s. 
"^Insufficient counts on target. 

''Elongated images due to tracking failure. 

■^Initial telescope focus not yet settled (3 points) and strongly aberrant ratio (1 point). 
^Prame transfer mode. 
^Comparison star saturated. 
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Table 4. Flux Values For New Transits of OGLE-TR-lllb " 



Mid-exposure (UTC) 


Mid-exposure (BJD) 


Flux 


Error 


2454574.520046 


2454574.52375 


0.9820188 


0.00153 


2454574.522041 


2454574.525745 


0.9797899 


0.00153 


2454574.522868 


2454574.526573 


0.9824249 


0.00153 


2454574.523734 


2454574.527439 


0.9802172 


0.00153 


2454574.524134 


2454574.527838 


0.9828871 


0.00153 



''Pull table available online. 



Table 5. Transit Parameters Based on Joint Fit to Eleven Light Curves 



Median value Formal Error^ Adopted Error 



Fitted Parameters 



k 


0.1261 


-1-0.0008, -0.0009 


+0.0010, 


-0, 


.0011 


a/Rt 


12.3 


0.2 


0.2 




i 


88.3 


+0.3, -0.2 


+0.3, 


-0, 


.2 




0.32 


0.03 


0.04 






0.252 


(fixed) 










0.30 


0.04 


0.05 






0.2582 


(fixed) 










0.6228 


(fixed) 








U2,V 


0.1587 


(fixed) 








Tc - 2453470 


0.56486 


+0.00065,-0.00067 


+0.00098, 


-0, 


.00089 


Tc - 2453787 


0.70934 


+0.00043,-0.00044 


+0.00131, 


-0, 


.00137 


Tc - 2453799 


0.75212 


0.00044 


+0.00117, 


-0, 


.00129 


Tc - 2454088 


0.79197 


0.00028 


+0.00071, 


-0, 


.00077 


Tc - 2454092 


0.80573 


0.00032 


+0.00080, 


-0, 


.00097 


Tc - 2454574 


0.54282 


+0.00040,-0.00039 


+0.00103, 


-0, 


.00068 


Tc - 2454578 


0.55469 


+0.00044,-0.00045 


+0.00061, 


-0, 


.00061 


Tc - 2454598 


0.62752 


0.00029 


+0.00055, 


-0, 


.00047 


Tc - 2454602 


0.64170 


0.00025 


+0.00044, 


-0, 


.00046 


Tc - 2454879 


0.63864 


0.00022 


+0.00044, 


-0, 


.00035 


Tc - 2454903 


0.72565 


0.00025 


+0.00029, 


-0, 


.00042 


Derived Parameters 












b 


0.35 


+0.04, -0.06 


+0.04, 


-0, 


.06 


Ti4 (sec) 


9647 


+47, -48 


+57, 


-60 


Rp {RjT 


1.019 


0.025 


0.026 




a (AU)'^ 


0.0473 


+0.0015, -0.0014 


+0.0015, 


-0.0014 



"Formal 68.3% credible interval from MCMC fit of all 11 light curves jointly. 

''Adopted error from residual permutation method, if greater than the formal error; 
see § 13.31 

'^Assuming R-t = 0.83 ± 0.02 Rq llSantos et al.ll2006h and using Rj = 71, 492fcm. 
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Table 6. Transit Parameters Based on Individual Fits to Each Light Curve 



Transit 








ri4" 














20050409 


0. 


127 ± 0.0047 


1.0 


9901 ± 381 




11.7 ± 1.2 


1.1 


88.1 ± 1, 


,5 


1.3 


20060221 


0. 


132 ±0.0038 


2.0 


9829 ± 170 


1.2 


12.2 ± 0.7 


1.4 


88.6 ±0, 


,8 


1.1 


20060305 


0. 


132 ±0.0051 


3.2 


9920 ± 234 


1.6 


11.8 ± 0.9 


1.6 


88.1 ±0, 


,9 


1.4 


20061219 


0. 


118 ± 0.0076 


3.1 


9436 ± 217 


1.0 


12.2 ± 0.8 




88.2 ±0, 


,9 




20061223 


0. 


120 ± 0.0064 


1.8 


9580 ± 255 


1.2 


12.2 ± 0.9 


1.1 


88.3 ±0, 


,9 




20080418 


0. 


129 ± 0.0036 


1.2 


8646 ± 1413^= 


3.2 


12.9 ± 1.4 


1.8 


87.7 ± 1, 


,5 


2.5 


20080422 


0. 


128 ± 0.0029 


1.2 


9837 ± 222 




12.9 ± 0.9 


1.2 


88.3 ±1, 


,0 


1.1 


20080512 


0. 


123 ± 0.0022 


1.9 


9482 ± 134 


1.5 


12.8 ± 0.6 


1.5 


89.0 ±0, 


,6 


1.0 


20080516 


0. 


127 ±0.0021 


1.5 


9677 ± 183 


1.4 


12.3 ± 0.7 


1.3 


88.6 ±0, 


,7 




20090217 


0. 


125 ±0.0025 


2.2 


9644 ± 193 


1.9 


12.4 ±0.7 


1.6 


88.7 ±0, 


,7 


1.1 


20090313 


0. 


127 ±0.0022 


1.9 


9390 ± 209 


2.7 


13.0 ±0.6 


2.2 


89.2 ±0, 


,6 


1.3 



^Formal individual MCMC fit value and error (scaled upward by factor / in adjacent column). 

''Factor by which the error in the previous column has been increased based on the residual 
permutation method; no value is given if the formal MCMC fit error was larger. See § 13.31 

'^The ill-constrained duration of the h alf-transit 20080413 was not used in any fits. 



Table 7. Timing Residuals 



Transit Number O-C (s) ct 



20050409=^ 


-155 


-271 ±83 


-3.3 


20050409*^ 


-155 


65 ±43 


1.5 


20060221 


-76 


8 ±75 


0.1 


20060305 


-73 


-41 ±82 


-0.5 


20061219 


-1 


-66 ± 103 


-0.6 


20061223 





-125 ± 114 


-1.1 


20080418 


120 


180 ± 87 


2.1 


20080422 


121 


-42 ± 52 


-0.8 


20080512 


126 


10 ±41 


0.3 


20080516 


127 


-13 ±46 


-0.3 


20090217 


196 


-1±37 


-0.02 


20090313 


202 


28 ±36 


0.8 



''Using our analysis of the original photome- 
try from Minniti et al. ( 2007). 

''Using the published time from rean alyzed 
photometry bv lPietrukowicz et al. 1 boidb . 



